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Abstract. A 2.29 GHz VLBI all-sky survey of ultra-compact radio sources has 
formed the basis of a number of cosmological investigations, which examine the 
relationship between angular-size and redshift. Here I use a sample of 468 such sources 
with 0.5 < z < 3.787, to investigate the isotropy of the Universe. The sample is divided 
into hemispherical sub-samples, over an all-sky 5° x 5° array, each of which is allowed 
to determine a value of f2 m , assuming that we are living in a spatially flat homogeneous 
isotropic ACDM model. If we regard the latter as a null hypothesis, then it fails the 
test - the results show significant anisotropy, the smallest value of f2 m being towards 
(7,6) = (253.9, 24.1)°, the largest in the opposite direction. This is close to the CMB 
dipole axis, but in the obverse sense. This is interpreted as meaning that the Universe 
is not spatially homogeneous on the largest scales, and is better represented at late 
times by a spherically symmetric model with a density enhancement at its centre. 
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1. Introduction 

In the large the Universe around us exhibits a remarkable degree of spherical symmetry, 
the most stringent constraint on anisotropy being provided by observations of the 
Cosmic Microwave Background (CMB). The near isotropy of the latter was apparent 
at the time of its discovery (Penzias & Wilson 1965; Wilson & Penzias 1967). Early 
COBE results showed a dipole component of AT = 3.335 mK (0.12 per cent), and 
a monopole temperature of 2.725 K (Kogut et al. 1993; Mather et al. 1999). After 
removing the dipole component the rms sky variation is AT = 45.3 fiK (0.0017 per 
cent) (see for example Percival et al. 2002) over the multipole range 2 < I < 1500. If 
we accept the Copernican principle then these observations mean that the Universe is 
homogeneous and isotropic; the corresponding Friedmann-Lemaitre-Robertson- Walker 
(FLRW) framework has been the the basis of most studies in observational cosmology. 
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However, while the Copernican principle remains untested, inhomogeneous models 
should not be dismissed (Clarkson & Maartens 2010; Ellis 2011). I report here a 
test of isotropy based upon the angular-size/redshift relationship, using ultra-compact 
radio sources as standard measuring rods; these objects have angular diameters in the 
milliarcsecond (mas) range, and linear sizes of order several parsecs. In fact the test 
reveals significant anisotropy, a tentative interpretation of which is that the Universe is 
not spatially homogeneous on the largest scales, and is better represented at late times 
by a spherically symmetric model with a density enhancement at its centre. Antoniou 
& Perivolaropoulos (2010) have already looked at Union2 Snla dataset in this context, 
which shows a similar anisotropy; my approach closely follows theirs. More recently, 
Longo (2012) has reported an anomaly in the angular distribution of quasar magnitudes, 
which is interpreted as evidence for a large-scale distant inhomogeneity. 



2. Ultra-compact radio sources 

The data to be used derive from an ancient VLBI survey of such sources at 2.29 GHz 
(Preston et al. 1985, hereafter referred to as P85). This survey employed a world- 
wide array of dishes, forming an interferometric system with an effective baseline of 
about 8 x 10 7 wavelengths; the survey gave total and correlated flux densities (fringe 
amplitudes) for 917 objects, S t and S c respectively. The ratio T = S c / S t is a measure 
of fringe visibility, from which angular size can be estimated (Thompson, Moran & 
Swenson 1986; Gurvits 1994): 

2 v / -lnrin2 

6= ZT> > I 1 ) 

7T-D 

where B is the interferometer baseline, in wavelengths. 

The potential of P85 in this context was first noted by Gurvits (1994), who 
considered a sample comprising 258 objects with redshifts z > 0.5. Using the same data 
set, Jackson & Dodgson (1997) extended this work to the full fl m /fl\ plane, finding 
best values fl m = 0.2 and Q\ = 0.8, if the Universe is spatially flat, later refined to 
Q m = 0.24, Q\ = 0.76 (Jackson 2004). The sample was updated with regard to redshift 
(Jackson & Jannetta 2006), to give 613 objects with 0.0035 < z < 3.787, of which 468 
have z > 0.5, which sub-set is used in this investigation. There are several reasons for 
ignoring sources with z < 0.5. As z falls below 0.5, the epoch of quasar formation comes 
to an end, and the nature of the population changes dramatically (see the cyan points 
in Fig. [TJ; there is a correlation between linear size and radio luminosity - the weaker 
sources are distinctly smaller, which would introduce an unacceptable selection bias. 
There is no evidence of such an effect when z > 0.5 (Jackson 2004); if it is there it is 
not pronounced. A residual bias should not appear to be anisotropic. 

A model of these objects, which gives an account of their status as standard 
measuring rods, is discussed in Jackson (2004). The model is supported by VLB A 
images (Kovalev et al. 2005), which show a bright compact 'core' at the end of a 
one-sided jet; the 'core' is believed to be the base of a continuous jet, rather than the 
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nucleus of the object. Beamed emission from the compact core dominates the observed 
structure. The importance of Doppler beaming cannot be over-stated; the rest-frame 
linear size of mas source components is known to be an increasing function of wavelength 
(Marscher & Shaffer 1980; Pearson & Readhead 1981); at first sight this would mean 
that the effective linear size should actually be a decreasing function of z. However, as z 
increases a larger Doppler boost factor D is required; it turns out that the ratio V/(l + z) 
is approximately fixed (Jackson 2004), so that the emitted frequency Vu T /(l + z) also 
fixed, where v x is the fixed reception frequency. Note that this ratio is not necessarily 
unity, but the fact that it is fixed means that the interferometric angular sizes upon 
which this work is based are fit for purpose. This is another reason for ignoring sources 
with z < 0.5, which appear to be non-relativistic. Again a residual bias should not 
appear to be anisotropic. 

A second point which needs clarification is resolution. With a baseline B = 8 x 10 7 , 
the Rayleigh resolution limit is about 2.6 mas. The mean angular size of the sources in 
question is 1.52 mas, somewhat below the Rayleigh limit. However, it is well-known that 
a simple interferometric technique can achieve a degree of super-resolution; the matter 
is discussed in detail in Kovalev et al. (2005). The limiting factor is signal-to-noise 
ratio SNR; for a simple Gaussian source the limiting angular size is 



spacing, thus b = (2B)^ 1 . As to the appropriate SNR, two figures are mentioned in 
Preston et al. (1985). The first is described as a systematic error of about 10 per cent; 
a systematic percentage error should not affect measures of fringe visibility, but I note 
that the corresponding resolution is 8n m = 0.20 mas. The second figure is an absolute 
random error of 0.2 Jy in the correlated flux density, in which case the SNR would vary 
from source to source; the median SNR value for those sources used in this invstigation 
is 17, giving 9\ im = 0.15 mas. These figures are indicative, but do confirm that lack of 
adequate angular resolution is unlikely to be an issue here. (However, for a radically 
different point-of-view see Pashchenko & Vitrishchak 2011). 

Fig. [l] is the angular-diameter/redshift diagram for the full sample of 468 sources 
discussed above. The continuous green line corresponds to the best-fitting spatially 
flat homogeneous isotropic ACDM model, characterized by two free parameters, the 
matter density Q m and the intrinsic linear size d. The best values are Q m = 0.29 and 
d = 7.76/i -1 pc (Hq = 100/i km sec -1 Mpc -1 ), being a fit to the unadorned unbinned 
data points; Fig. [2] shows the corresponding confidence regions (the green curves). 

I have undertaken a similar analysis of the Caltech-Jodrell Bank 5 GHz survey, 
giving similar results, which survey cannot however be used in this context because it 
is not an all-sky one (Taylor et al. 1996; Jackson 2008). 




(2) 



where b is the half-power beam width, which I have taken to be be half the fringe 
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Figure 1. Angular-diameter/redshift diagram for the full sample of 468 sources (blue 
points), and the corresponding best-fitting curve f2 m = 0.29 and d = 7.76/i _1 pc. The 
blue dashdot line is an estimate of the resolution limit #n m = 0.15 mas. The 145 cyan 
points have z < 0.5, and are not used in this investigation. 

3. Anisotropy 

To test isotropy, I have considered hemispheres, typically containing 234 ±11 sources, 
which data points are allowed to determine fl m and d; the values so determined now 
depend upon the particular hemisphere, defined by the Right ascension and declination 
(a, S) of its axis. The coordinates used here are B1950, as listed in P85. Fig. [3] is a 
pseudocolor Mercator projection based upon an evaluation of Q m on a 5° x 5° grid over 
the range 0° < a < 360°, —60° < 5 < 60°. The plot shows a pronounced asymmetry, a 
measure of which is the parameter D = AQ m evaluated over opposing hemispheres. The 
global maximum value is D max = 0.626 - 0.135 = 0.491 towards (I, b) = (253.9, 24.1)°, 
which corresponds to the smaller value of fl m , indicated by the cyan V mark in Fig. [3j 
the opposite direction is indicated by the grey V mark. 

The red and blue x marks indicate the CMB dipole in the Local Group frame, with 
velocity 627 ± 22 km sec" 1 , the 'hot' direction being (I, b) = (276 ± 3, 30 ± 3)° (Kogut 
et al. 1993). The proximity of these two directions is quite striking; however, the CMB 
dipole cannot be attributed to a peculiar velocity induced by gravitational attraction 
towards the region of enhanced density, which would be in the wrong direction, a 
point to which I shall return later. The features indicated by the A and square 
marks are for later reference. The figure reported by Antoniou & Perivolaropoulos 
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Figure 2. Confidence regions, one and two a: green corresponds to the the full 
sample, blue and red to opposing extreme hemispheres. These are marginalized, so 
that projection on to each axis gives confidence intervals for the respective parameter. 

(2010) is L> max = 0.30 - 0.19 = 0.11 towards (I, b) = (309+jg, 18^°)°, which again 
corresponds to the smaller value of Q m . The Union2 Snla dataset covers the redshift 
range 0.015 < z < 1.4 (Amanullah et al. 2010). 

The blue and red confidence regions in Fig. [2] correspond to the hemispherical 
samples which determine (fi m )min and (fi m ) max , which samples I shall call S min and 
S max . I have considered one way in which the differences shown in Fig. [2] might be 
spurious. The fitting problem is close to degeneracy with respect to the parameters 
Q m and d, because we cannot use nearby sources, which would otherwise fix the latter 
(Jackson & Dodgson 1996). I have considered the possibility that there are in fact no 
significant differences between S min and S max , other than small accidental ones which 
due to the degeneracy are nevertheless large enough to bring about the divergences seen 
in Fig. |2j This possibility is discounted by the following considerations; I have further 
divided the two hemispherical samples, into S min (z < 1.5) and S min (z > 1.5), and similarly 
for S max . It is indeed the case that the low-redshift sub-samples are virtually identical, 
but in the high-redshift case there is a statistically significant difference: the two mean 
angular-sizes are 9 m i n (z > 1.5) = 1.35 ± 0.05 mas and 9 max (z > 1.5) = 1.56 ± 0.06 mas 
(1 cr errors), which difference is not accounted for by a difference in the respective values 
of z. The two values of 9(z > 1.5) are the driving force behind the results presented 
here. 

I have considered sensitivity to the exclusion of high-redshift sources, by looking at 
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Right ascension 

Figure 3. Pseudocolour plot showing the distribution of fl m over the sky; full sample 
of 468 sources, 0.5 < z < 3.787. There are no prominent features in the polar regions. 
Crosses indicate the CMB dipolc in the Local Group frame; the other symbols delineate 
prominent features for further reference. 

restricted samples defined by 0.5 < z < z max - The pattern remains reasonably stable as 
long as 2 max > 2.0 (see Fig. [5]); at z max = 1.5 there is little trace of the original structure. 
The results presented here thus do not depend upon a few high-redshift objects. 

Could statistical fluctuations in the distribution of for example high-redshift 
sources between opposing hemisperes account for the apparent asymmetry? I have 
evaluated D max over an ensemble of mock samples, generated by randomizing source 
positions while leaving the corresponding redshift and angular-size data unchanged. 
The distribution of -D max values so produced is broad, with a median value of 0.508. 
The above value -D max = 0.491 is thus inconclusive. As an alternative test, I have 
first subjected the full list of 468 sources to a random shuffle, and then divided the 
list into two independent sub-samples, the first comprising all even-numbered members, 
the second all odd-numbered such members; Figs. [5] and [6] show the hemispherical 
distribution of Jl m over the sky for the even and odd sub-samples respectively. If the 
prominent features delineated in Fig. [3] are due entirely to statistical fluctuations, then 
Figs. [5] and [6] would not be expected to show the same features, whereas it is quite 
clear that they do. Figures [5] and [6] are very similar; they are highly correlated, with 
a correlation coefficient of 0.413. Remembering that the even and odd sub-samples 
have no members in common, the probability of this situation arising by chance is less 
than 1 x 10~ 9 (see for example Freund & Walpole 1980). I conclude that the observed 
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Right ascension 

Figure 4. Pseudocolour plot showing the distribution of f2 m over the sky; restricted 
sub-sample of 380 sources, 0.5 < z < 2.0. The features highlighted in Fig. [3] have not 
been eliminated by dropping the high-redshift sources. 

asymmetry is intrinsic. There are 468! distinct permutations of the original list, and I 
have examined a modest number of these; Figs. [5] and [6] represent a typical example. 

The remaining possibility is that the apparent anisotropy is an instrumental effect, 
affecting all of the above-mentioned samples and sub-samples in like fashion. In this 
respect a suspect feature of the full sample of 468 sources is that not all of the flux 
measures derive from the same system; 339 are pure P85, in that they have both 
correlated and integrated flux densities listed in P85; the remaining 129 have correlated 
fluxes in P85, without the corresponding integrated fluxes; in these cases the latter 
were taken from the Parkes catalogue PKSCAT90 (Wright & Otrupek 1990). As the 
correlated fluxes are the critical ones, I do not regard this as a serious deficiency; to 
confirm this I have discarded the 'mixed' sources and repeated the primary computation 
described above: Fig. [7] shows the Q m distribution determined by the set of pure P85 
sources, which is very similar to Fig. [3} The anisotropy is clearly not an artefact of 
instrumental inhomogeneity. Finally there is the matter of uniformity of sky cover. 
Without the additional Parkes sources the pure P85 sample is deficient in southern- 
hemisphere ones (5 > —41.9°); with the additional Parkes sources the sample is uniform, 
the mean number per hemisphere being 233.95 ± 11.16, very close to the binomial figure 
234 ± 10.82. 
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Right ascension 

Figure 5. Pscudocolour plot showing the distribution of O m over the sky; even sub- 
sample of 234 sources, see text. 
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Figure 6. Pseudocolour plot showing the distribution of Q m over the sky; odd sub- 
sample of 234 sources, see text. 
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Right ascension 

Figure 7. Pseudocolour plot showing the distribution of fl m over the sky; restricted 
sample comprising 339 pure P85 sources, 0.5 < z < 3.787. 

4. Interpretation and a toy model 

We cannot of course maintain that the Universe is well-represented by a flat 
homogeneous isotropic ACDM model, but that a different such model is required 
according to direction. Nevertheless, I believe that these results mean that there is 
more dark matter in some directions than in others, whatever the geometrical setting. 
To see this, we consider a source represented by a small plane circular disc, perpendicular 
to the light-ray going from its centre to the observer. In vacuum, the bundle of rays 
going from the edge of the disc to the observer defines a cone with the observer at its 
apex; in the non-vacuum case the cone is refocussed by the gravitational attraction of 
the matter within it, as we trace the rays back towards the source, which refocussing has 
a magnifying effect. Indeed if the source is sufficiently distant then the cone begins to 
reconverge, and the apparent size of the source begins to increase (see for example Ellis & 
Tivon 1985). A reasonable explanation of the fact that the sources in the hemispherical 
sample S max appear to be systematically larger than those in S min , is that the rays in the 
direction of the former are passing through denser matter than those in the direction of 
the latter. In one or two simple cases an exact analytical treatment is possible (Jackson 
& Dodgson 1996), based upon the Ehlers-Sachs equation (Ehlers & Sachs 1959, Sachs 
1961, Pirani 1965), which rederives some of the elementary cosmological results relating 
to distance. 

The matter concentration revealed here cannot be behaving like a simple attractor, 
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because the motion of the Local Group would be towards the latter, whereas it appears 
to be in the opposite direction. However, this is exactly what an off-centre fundamental 
observer (at rest relative to the local Hubble flow) would expect to see in a spherically- 
symmetric dust-filled inhomogeneous model with a central concentration (Raine & 
Thomas 1981; Maartens et al. 1996a; Maartens et al. 1996b; Humpreys, Maartens 
& Matravers 1997). According to this picture the apparent motion of the Local Group 
is due largely to a peculiar cosmological redshift, rather than a peculiar flow (Paczynski 
& Piran 1990). 

Large-scale coherent peculiar motions would be a feature of such a model, which 
must be compared with observations. The local such flow is believed to be due to the 
Great Attractor, a massive concentration of elliptical galaxies apparently at rest in the 
CMB frame, at a distance of 43.5/i" 1 Mpc towards (I, b) = (307, 9)° (Lynden-Bell et al. 
1988); this result was based upon a sample of 400 elliptical galaxies closer than 80/iT 1 
Mpc. Infall (rather than uniform bulk flow) gives a good account of the near-side flow; 
however, evidence for infall from the far side of the attractor is inconclusive (Dressier 
& Faber 1990). Watkins, Feldman & Hudson (2009) and Feldman, Watkins & Hudson 
(2010) examined a compilation of similar samples, and concluded that the mean peculiar 
motion of galaxies within a spherical volume of radius c. 100/i _1 Mpc is 416 ± 78 km 
s -1 , towards (/, b) = (282 ±11, 6 ± 6)°, and that this velocity appears to be increasing 
with distance. As the authors remark, a flow of this amplitude on such a large scale is 
not expected in the WMAP5-normalized ACDM cosmology. 

Kashlinsky et al. (2008, 2009, 2010) have looked for CMB temperature fluctuations 
induced by the kinematic Sunyaev-Zel'dovich effect, from which the velocities of the 
associated clusters can be deduced. These authors report a positive detection with bulk 
velocities in the range 600 to 1000 km sec -1 towards (l,b) = (283 ± 14, 11 ± 14)°, with 
most of the signal arising from a shell between 150/i _1 and 600/i~ 1 Mpc; such flows are 
not easily accomodated within the standard inflationary ACDM model (Atrio-Barandela 
et al. 2010). However, the latter results have been discounted by others (Keisler 2009; 
Osborne et al. 2011; Mody & Hajian 2012; Hand et al. 2012). 

For purposes of comparison, I have considered a simple spherically symmetric zero- 
energy Newtonian model: 

2 2GM 1 

HM) + 3 Ar(M) ' (3) 

where r(M) is the comoving radius of a sphere of mass M centred on the point 
of maximum concentration, and v(r) is the radial velocity relative to the centre. 
Introducing p as the mean density within a sphere of radius r, equation ^ becomes 

2 MGpr 2 1 A 2 , A , 
v 2 = 3 H + -Ar 2 . (4) 

We consider a snapshot and radial power series expansion about a point at r^. 

v = v + (r-r )(^) +l(r-r ) 2 (^] + (5) 



dr / n 2 V dr 



o 
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The first-derivative in equation ^ gives the local Hubble flow: 
v H = H Ar = [-AnGA- 1/2 (p - p) + A 1 ' 2 ] Q Ar, 



(6) 



where Ar = r — r and A = 8irGp/3 + A/3. The second derivative gives a local radial 
peculiar flow; on a shell of radius Ar, this is (dropping the subscript): 



.4 



-l 

— [4vrG(p-p)] 2 + A~ 1/2 4ttG 



- (P-P) +"T 
r ar 



Ar 2 



(7) 



A positive value for v B corresponds to peculiar motion away fron the centre. We are 
considering distributions in which p—p > and dp/dr < 0, so that either sign is possible; 
for example, if p — p » | dp/dr | and p and p are not too different, the dominant term 



is 



A- 1/2 47rGp 



p — p\ Ar 



P 



p — p\ Ar 
P 



-H Ar. 



(8) 



Taking (p — p)/p and Ar/r to be 0.20, Ar = 100 Mpc and H = 70 km sec -1 Mpc -1 , 
equation ^ gives v B = 140 km sec -1 . The main point of this rough estimate is to 
show that the model is not obviously at variance with the developing observational 
evidence, and that equations ^ and ^ are sufficiently flexible to accomodate possible 
contingencies. 
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